/ 



THE UNDERWATER APPLICATION OF 
EXOTHERMIC WELDING 



Arthur H. Anderssen 



f 



1 



THE UNDERWATER APPLICATION OF 
EXOTHERMIC WELDING 
by 

ARTHUR H. ANDERSSEN 

// 

B.S., Auburn University 

( 1962 ) 

SUBMITTED IN PARTIAL FULFILLMENT 
OF THE REQUIREMENTS FOR THE 
DEGREE OF OCEAN ENGINEER 
at the 

MASSACHUSETTS INSTITUTE OF 
TECHNOLOGY 



June, 1972 



M.I.T. THESES ARE THE PERMANENT PROPERTY OP THE 
INSTITUTE AND SHALL NOT BE REPRODUCED EITHER IN 



PULL OR IN PART. 



2 



The Underwater Application of 
Exothermic Welding 

by 

Arthur H. Anderssen 

Submitted to the Department of Ocean 
Engineering on May 12, 1972, in partial 
fulfillment of the requirements for the 
degree of Ocean Engineer. 

ABSTRACT 

Existing welding processes which have been adapted 
for underwater use impose serious problems and often produce 
welds of questionable quality. This thesis presents an in- 
vestigation of a process previously untested in the under- 
water environment. The overall objective of the study is 
the conceptual design of a device which could be used to 
exothermicly weld a stud, bar, or padeye to underwater 
structures. The process is commonly called thermit welding. 

Several problem areas of both long and short range 
interest are considered. A theoretical temperature simula- 
tion model is developed and verified. The effects of wetness 
on weld strength and the ability of the thermit process to 
weld on a "dirty" surface are investigated experimentally. 

A conceptual design is prepared using experimental results. 

By far the most important conclusion arising out 
of the study is that the thermit process can be applied 
successfully underwater. 
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DEFINITION OF SYMBOLS USED 



Definition 

Specific heat of saturated 
liquid 

Constant dependent upon 
surface fluid combination 

2.71828 

Conversion factor 

Gravitational accelleration 

Convective heat transfer 
coefficient 

Latent heat of vaporization 

Bessel Functions 

Thermal conductivity 

Length 

Heat flux 

Radius coordinate 

Time 

Temperature 

* 

Depth coordinate 
Thermal difusivity 
Coefficient of expansion 
Angle 

Dynamic viscosity 
Density 

Surface tension 



Units 

Btu/lbm-°F 



lbm-f t/lbf-hr 2 
ft/hr 2 

Btu/hr-ft 2 -°F 
Btu/lbm 

Btu/hr-f t-°F 
in or ft 
Btu/hr-ft 2 
in 

sec or hr 
°F or °R 
in 

ft/hr 

1/°F 

degrees or radians 

Ibm./ft-sec 

lbm/ft 

lbf/ft 
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I . INTRODUCTION 

A. THE OVERALL PROBLEM AND ITS SETTING: 

Man's capability to work in the ocean environment 
has increased dramatically. This increased capability has 
created a new need for cutting and joining techniques which 
can be applied in the underwater environment. Thus, we see 
an increased amount of goal directed research around the 
world in underwater cutting and joining techniques. For 
example, the petroleum industry's offshore wells have 
created a need for suitable joining procedures. This need 
has contributed to the development of underwater hyperbaric 
welding techniques and hardware over the past few years. 

However, underwater arc and oxyacetylene welding 
using a diver, today and for the foreseeable future, appears 
limited to depths of less than 1000 feet, even employing 
recent developments in saturated diving. Yet requirer ~r.t s 
exist and more will soon develop for welding in depths far 
greater than this. Considering the fact that about 95 per- 
cent of the entire ocean floor is deeper than 1000 feet, 
there is a strong need for developing a capability for 
welding in the deep sea, particularly for salvage work. One 
way to perform the welding work is to use a deep submersible 
and employ a welding process which can be controled remotely 
The use of an exothermic welding process (thermit) may offer 
a possible solution. It is a relatively simple process 
requiring no power supplies, and could be remotely 
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accomplished once the welding device was positioned. 

The possible value of such a device has been 
recognized particularly by U.S. Navy personnel engaged in 
salvage work. As a result, toward the end of 1970, the 
Navy 'granted funds to the Massachusetts Institute of 
Technology to investigate characteristics of the thermit 
reaction in the underwater environment and its applicability 
to the remote-controlled welding of a 50-ton padeye to 
underwater structures. The results of this research 
essentially led to the conclusion that the application of 
exothermic welding processes in the underwater environment 
deserved further study. This thesis undertakes that study 
with particular reference to a means of welding a padeye or 
stud to a deeply submerged object for salvage purposes. 

Before stating more specific objectives, however, 
the reader might find a brief introduction to thermit welding 
helpful . 
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B. EXOTHERMIC WELDING: 

Although its use is not as universal as arc 
welding, for example, thermit welding has been used with 
considerable success in some applications above water. 
Consequently, a body of knowledge has evolved which must be 
understood by those attempting to apply the process to a new 
environment (such as underwater). For this reason, a brief 
summary of thermit welding processes and major properties 
of resulting welds is included here. 

When metal oxides having low heats of formation 
are reacted with reducing agents which, when oxidized, have 
high heats of formation, the reaction is exothermic. The 
terminology "thermit welding" refers to the group of welding 
processes which use the superheated metal and slag produced 
by the above reaction to achieve coalescense with or without 
the application of pressure. The following are typical 
thermit reactions: 



Reaction 




» 














Theoretical Temp. 
Achieved 


3 Fe 3 0 4 


+ 


SA 1 


-> 


SFe 


+ 


4A1 2 C 3 


+ 


719.3 


Heal 


5590 c F/3088°C 


3FeO 


+ 


2A1 


-> 


3Fe 


+ 


A1 2 0 3 


+ 


187.1 


Kcal 


4532 o F/2500°C 


Fe 2 0 3 


+ 


2A1 


-* 


2Fe 


+ 


A1 2 0 3 


+ 


181.5 


Kcal 


5360 o F/2960°C 


3CuO 


+ 


2A1 




3Cu 


+ 


ai 2 o 3 


+ 


275.3 


Kcal 


8790°F/il865 o C 


3Cu 2 0 


+ 


2A1 


-> 


6 Cu 


+ 


A1 2 0 3 


+ 


260.3 


Kcal 


5680°F/3138°C 


3NiO 


+ 


2A1 


-y 


3Ni 


+ 


A1 2 0 3 


+ 


206.6 


Kcal 


57^0°F/3171 o C 


Cr 2 0 3 


+ 


2A1 




2Cr 


+ 


A1 2 0 3 


+ 


5-46.5 


Kcal 


5390°F/2977°C 
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Reaction 
3MnO + 2A1 
3Mn0 2 + 4A1 



3Mn + A1 2 0 3 + 403.0 Kcal 

3Mn + 2A1 2 0 3 + 1041.0 kcal 



Theoretical Temp 
Achieved 

4400°F/2427°C 

5020°F/2771°C 



The first reaction above is by far the most com- 
mon, and the major portion of my discussion will refer to 
it. The temperature required to initiate the reaction is on 
the order of 2200°F/1205°C . This temperature is generally 
provided through an ignition powder which, in turn, is 
ignited by a spark, electric match, etc. Although the 
theoretical temperature resulting from the reaction is 
5590°F/3088°C , radiant heat losses and losses to the reac- 
tion vessel reduce this temperature to about 4600°F/2538°C . 
In practice, other additions and impurities reduce the 
temperature of the filler weld metal to about 3800°F/2093°C . 

There are several methods of employing thermit in 
joining. The most common is fusion welding where the molten 

metal produced by the thermit reaction is allowed to flow 

« 

between the parts to be joined. The thermit is reacted in 
a crucible and is tapped into a mold after the reaction is 
complete and slag separation has occurred. A portion of the 
molten metal is contained between the parts by the mold and 
a portion allowed to flow into an overflow chamber or riser. 
Flow through the mold performs two important functions: 
surface cleaning and preheating. Welds of various sizes and 
configurations can be accomplished in this manner with a 



12 



variety of mold types. Figures 1 through 3 illustrate 
various hardware for different applications. 

Pressure thermit welding refers to the process 
where the molten metal and slag are used to provide suffi- 
cient heat to join the parts themselves without adding 
material to the weld. In this method, the thermit compo- 
sition is controlled so that slag with a high freezing point 
is produced. The thermit is reacted in a separate crucible, 
configured in such a way that slag is tapped first into a 
volume surrounding tightly fitted parts. A frozen layer of 
slag is then formed next to the parts to be joined before 
the molten metal enters the mold. The heat content of the 
molten metal and frozen slag brings the temperature of the 
parts up to forging level, and sufficient pressure is 
applied to produce a bond. 

A third, more recently applied method of using 
thermit for joining is thermit brazing. As in pressure 
thermit welding, thermit mixtures are used solely for their 
heat content. Brazing fluxes and filler metal are positioned 
so that when the parts to be joined are exothermicly heated, 
the flux will clean the surfaces and the brazing material 
will flow between the parts by capillary action. 

The properties of thermit welds are affected 
mainly by thermit composition and cooling rate. Unfortun- 
ately, since most research has been performed by companies 
having a vested interest in the results, little hard data is 
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Figure 1 



THERMIT WELDING KIT FOR JOINING RAILS 




.. 



14 




Figure 2: TYPICAL GRAPHITE MOLD FOR JOINING COPPER CABLE 
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Figure 3: HALF SECTION OF BAR TO PLATE THERMIT WELDING KIT 
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available. Composition of the weld metal is controlled by 
adding metallic elements either in the form of pieces or 
compounds which produce secondary exothermic reactions. The 
A.W.S. Welding Handbook reports that as cast tensile 
strengths up to 130,000 psi, and elongations from almost 
zero to 49 percent are possible by controlling composition. 
Thermit weld strengths dr. /’airly easily matched with 

the material welded. figures m and 5 illustrate this point 
through the results of strength tests performed on welded 
and unwelded reinforcing bars by Thompson and Lichtner Co., 
Inc . 



If the thermit is free of alloying elements (i.e., 
only Pe 2 0 3 and A1 are used), the effect of aluminium on weld 
properties can be determined. As would be expected, excess 
aluminium serves to reduce porosity, but contributes to 
brittleness. Silicon has the same effect and is a major 
cause of brittleness whenever sand molds are used. Manganese 
also promotes soundness in amounts of about 1 percent and 



does not have the severe effect cn brittleness as dc 

aluminium and silicon. 

The cooling rate has a significant influence on 
resulting hardness. In past research, the exact nature of 
which is proprietary, hardness of the fusion zone correlated 
positively with cooling rate. Figure 6 illustrates this 
effect although exact cooling rate was not determined. 



STRESS IN 
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Figure 4: RESULTS OF TENSION TEST ON 18S- INTERMEDIATE 



GRADE REBARS, SOURCE: THERMEX METALLURGICAL, INC . 



STRESS IN 



18 




Figure 5: RESULTS OF COMPRESSION TEST ON 18S-INTERMEDIATE 

GRADE REBARS, SOURCE: THERMEX METALLURGICAL, INC. 



HARDNESS: ROCKWELL "B 
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Figure 6: EFFECT OF COOLING RATE ON HARDNESS 



FASTER COOL MEDIUM COOL SLOWER COOL 
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C. OBJECTIVES AND APPROACH: 

The overall objective of this thesis is the 
conceptual design of a device which could be used to exo- 
thermicly weld a stud, bar, or padeye to underwater struc- 
tures. As an obstacle to reaching this objective, the 
underwater environment contributes to five specific problem 
areas not commonly ■e-red in above-water exothermic 

welding . 

First, the chemical composition of the environment 
is vastly different from that encountered on the surface. 

The presence of water -or moisture in the vicinity of the 
joint to bp welded or the reaction chamber, will affect the 
properties of the weld and the characteristics of the reac- 
tion. Wet welds made using other welding processes adapted 
for underwater use have been shown to be subject to extreme 
stress corrosion cracking and hydrogen embrittlement as a 
result of the environment. The objective related to this 
problem is to determine just how much water or moisture can 

f 

be tolerated between the surfaces to be welded ar.d still 

produce a sound weld. The approach to this problem is 
experimental in nature. It is assumed that the thermit 
reaction itself must take place in a completely dry environ- 
ment . 

Pressure is a second problem area encountered. 
Although pressure is one reason that a remotely controllable 
process such as exothermic welding may be beneficial for 
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underwater application, device configuration will certainly 
be influenced by the factor. The objective related to this 
problem is to determine a means of providing a proper pres- 
sure balance between the reaction vessel and the environment 
at the time of tapping the molten metal. The approach to 
this problem is analytical in nature. 

Surface contamination, although somewhat of a 
problem above water, is a third problem area underwater. 
Contamination may vary from minor surface corrosion to 
painted surfaces to heavy marine growth. Obviously, at the 
time of welding, the surfaces to be joined must be free of 
foreign material, implying the need for pre-cleaning by the 
thermit welding device itself or a separate operation. The 
objective related to this problem is to determine that 
degree of surface contamination which can be effectively 
removed by flow of the molten metal itself between the sur- 
faces to be joined. The approach to this problem is 
experimental in nature. 

A fourth problem area is that of post it le water 

flow velocities in the area where the weld is to be accom- 
plished. Again, device configuration is affected since some 
means of holding the device in position while welding is 
accomplished must be devised. However, this problem is 
believed to be minimal in comparison to the others and is 
only treated conceptually during the design phase of this 
study . 
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Finally, the environment itself forms a large heat 
sink, affecting cooling rates and ultimately, the properties 
of the weld. Cooling rates and peak temperatures are of 
particular interest if thermit welding is found to be a 
viable process for underwater use, since metallurgical 
changes and weld properties can be predicted with knowledge 
of temperature history. Such information might also be 
useful in design stages if cooling rates could be altered by 
device configuration. Since the cooling rate problem could 
have both short and long term impact, a large portion of the 
thesis is devoted to it. The objective related to this 
problem area is to develop and verify a model of the thermit 
bar-to-plate welding process such that the temperature 
history of the metal could be predicted. The approach is 
analytical with experimental verification. 
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D . SUMMARY : 

Chapter II covers the development of the temper- 
ature simulation model. The exact solution to the heat flow 
equation is obtained, but since boundary conditions do not 
permit the evaluation of integration constants, a model 
using finite difference techniques is constructed. 

The experimental procedures are discussed in 
Chapter III. Data involving temperature histories, wetness 
level, and cleaning action was obtained during the experi- 
ments. Then, theoretical and experimental results are 
presented and compared in Chapter IV. 

Results are utilized in Chapter V to develop a 
conceptual design of an underwater welding device. 

Finally, conclusions and recommendations for 
further research are presented in Chapter VI. 
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II. THEORETICAL MODEL FOR HEAT TRANSFER 
A. OVERVIEW: 

Unlike most welding processes, the heat source in 
thermit welding has a well defined shape since a mold is 
used which surrounds the parts to be welded. In addition, 
the device itself does not move during the welding operation. 
These two features of thermit welding allow the researcher 
great flexibility when constructing a model of heat trans- 
fer: a coordinate system can be chosen which facilitates 

description of the heat source, and coordinate movement need 
not be provided for. 

A need to weld bars on thick plates (greater than 
1/8 inch) is expected. Therefore, plate thickness must be 
included if any accuracy in predicting temperature history 
is to be expected, implying the need for a three dimensional 
model if Cartesian coordinates were selected. However, 
since this thesis is concerned with welding a bar or stud to 
a plate, the symmetry cf the problem permits the selection 
of a cylindrical coordinate system. Only the dimensions !, r" 
(radius) and "z" (thickness) can be used and still provide 
a complete description of heat flow. 

Temperatures are expected to vary from approxi- 
mately 40°F to 4000°F. For this reason, any model must 
consider material properties variable, specifically, thermal 
conductivity (k) and thermal dif'fusivity (a) in this case. 
Furthermore, in the model developed here, heat losses during 
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solidification are also taken into account. Since the 
welding is to be performed underwater, heat losses attribu- 
table to the environment must be considered. These include 
convection, boiling and radiation. 

The following sections of this chapter present the 
governing differential equation, the finite difference re- 
duction, the boundarv conditions necessary for solution, and 
finallv, the overall computer model. 
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B. HEAT FLOW EQUATION: 



The general equation of heat flow without source 



terms and assuming conductivity to be independent of position 

. 1 
is : 



V 2 T 



1. 9T 
a 9 1 



where a = — 
cp 



( 2 . 1 ) 



In cylindrical coordinates, if T is independent of 0, the 
equation appears as follows in partial differential notation: 



9 2 T 1.9T 9 2 T = l.JT 

9r 2 r 9r 9z z a t 



( 2 . 2 ) 



By separating the variables, the following solution can be 
2 

obtained : 



T = [ A • e Y at ]-[B.J 0 Ur) + OYoUr)] 

r -p, -iz/y 2 -c 2 , u +iz/y T -C 7 i 

• [D‘e + E*e J 



(2.3) 



where A , P 3 c } n , p 5 y end £ are an v s v:hieh rpnet he 

evaluated by applying the appropriate boundary and initial 
conditions . 

Unfortunately, the boundary conditions in this 
case are such that the constants cannot be evaluated. 
Therefore, the equation cannot be solved exactly, and a 
different approach must be used to predict temperature his- 
tories . 



"hee page 7 for definition of symbols. 

2 

See Appendix A for detailed solution procedure. 
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C. FINITE DIFFERENCE SOLUTION: 

Finite difference techniques can be applied to an 
unsteady-state conduction problem as long as the initial 

q 

temperature distribution and boundary conditions are known. 
This method is especially convenient when a digital computer 
is available. 

First, a grid which describes the body in question 

is laid out as in Figure 7* A finite difference equation 

is then derived from the differential equation and is as 
!| 

follows : 

k* (T] - T 0 )«2‘TT«(r - r/2)*Az k • (T 2 - Tp ) • 2 • tt • r • Ar 

Ar Az 

k*(T 3 - To) , 2 , iT , (r + r/2)*Az , k • (T 4 - T p ) • 2 • tt • r • Ar 

Ar Az 

p • c • 2 • tt *r • Ar • Az • (T 0 - To) ,, \ 

At 

The above equation holds for any interior point in a two 
dimensional grid described by cylindrical coordinates (r 
and z). The equation is solved for the temperature (Tq) at 

1 

one time interval (At) later and is as follows if Ar = Az = s 



"5 

See Kreith, Frank, Principles of Heat Transfer, 
International Textbook Co., Scranton, Pennsylvania^ 1969; 
or Rohsenow, Warren M., and Choi, Harry, Heat , Mass , and 
Momentum Transfer, Prentice-Hall, Inc., Englewood Cliffs, 
New Jersey - ^ 1965 for discussions of finite difference 
methods . 



For clarity, the equations in this section are 
written using numerical subscripts rather than alphabetic 
as in the computer program. 
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Figure ?: CYLINDRICAL GRID FOR FINITE DIFFERENCE COMPUTATION 
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To = i-[Ti-(l - s/ ( 2 • r ) ) + T 2 + T 3 *(l + s/(2-r)) + T 4 ]+ 
(1 - 4/M ) -To 



where M = 



a • At 



In this case, M must be greater than 4 for the solution to 
converge. This infers that the selection of At is dependent 
upon grid size and the material used. 

By applying the equation to every nodal point in 
the network, the temperature distribution at time t + At can 
be determined from a given distribution at time t. Repeated 
application of the equation produces the temperature distri- 
bution at time t + n*At, where n is an integer. 

The equations for points located on the boundary 
of the body are somewhat different and depend upon the 
boundary conditions, i.e., heat losses. These are discussed 



in the next section. 
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D. HEAT LOSSES: 

In this model, heat losses at the plate boundary 
are caused by convection, boiling, and radiation. However, 
unlike underwater arc and oxyacetylene welding, in this 
case, convection accounts for the majority of heat transfer 
to the environment. This is due to the thermit mold forming 
an effective insulation between areas of extremely high 
temperature and the environment. 

When and if the surface temperature of portions 
of the plate in contact with water exceeds saturation tem- 
perature, nucleate boiling is also possible. This condition 
exists until the surface temperature of the plate reaches a 
point where an unstable vapor film covers the surface. At 
this point heat losses due to nucleate boiling reach a 
maximum; the point is sometimes referred to as that at 
which peak heat flux occurs. Subcooling of the environment 
tends to increase both the temperature of this point and the 
peak heat flux value, although literature search produced no 
studies which’ were directly applicable to this problem (pool 
boiling in subcooled water). However, the temperatures at 
the plate/water interface were not expected to reach peak 
heat flux level, and therefore, this point and further 
regimes of boiling (transition and film were not considered 
in this model. 

When dealing with heat losses from a flat plate 
in the underwater environment, the orientation of the plate 
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is particularly important. As in most previous research on 

heat transfer in underwater welding, this study considers 

5 

only horizontal plates. But even in this restricted case; 
heat loss mechanisms on the bottom surface are different 
from those on the top. Differences in convection phenomena 
are clear. Heat fluxes attributable to convection are 

g 

derived in Appendix B and summarized as follows: 

Top Plate Convection: 

q c /A = 32.2(T plate - T wat0r )‘ /3 [BTU/hr/ft 2 ] (2.6) 
Bottom Plate Convection: 

q c /A = 12.2(T p i ate - T„a ter .) S/ ‘ [BTU/hr/ft 2 ] (2.7) 

For nucleate boiling on the top surface, the 
correlation of experimental data presented by Rohsenow and 
Choi was considered to be valid within the accuracy required 

7 

by this model. Heat flux on the plate's top surface attri- 
butable to nucleate boiling is derived in Appendix B and is 
as follows: 

Top Plate Nucleate Boiling: 

q b /A = 1.7494(T plate - T sat )’ [BTU/hr/ft 2 ] (2.8) 

5 

For example, see Staub, James A., Fr., Temperature 
Distribution in Thin Plates Welded Underwater , Unpublished 
M.I.T. Engineers Thesis, June, 1971* 

g 

°See Kreith, og_. cit . , pp. 

7 

'See Rohsenow and Choi, o£. cit . , pp . 211 - 237- 
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However, this study hypothesizes that the above correlation 



horizontal, bubble activity on the bottom surface is restric- 
ted, i.e., the normal boiling mechanism of bubble formation, 
collapse and ascent does not occur. What is hypothesized 
to occur is the formation of a semistable vapor film which 
effectively restricts heat transfer to a radiation mode. 

Since temperatures were not expected to be high enough for 
radiation to have any significant effect, the phenomenon is 
modeled as zero heat loss from areas on the bottom surface 
which are at a high enough temperature (T , ) to produce 

S 3.0 

the vapor film. In equation form: 

Bottom Plate: 



. g 

does not hold for the bottom surface. Since the plate is 



q b /A = 0 
q r /A = 0 



( 2 . 9 ) 



( 2 . 10 ) 



The total heat loss from either surface is simply 



the sum of the individual components as follows: 




( 2 . 11 ) 



Based in part on discussion with Warren M Rohsenow. 
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E. OVERALL MODEL AND COMPUTER SIMULATION: 

As stated previously, since it is envisioned that 
a thermit welding device designed for underwater stud welding 
would be cylindrical in shape, at least in the area where 
it contacted the plate, cylindrical coordinates are used in 
the model. The thermit mixture is assumed to be reacted in 
a separate crucible and tapped into a space bounded by the 
plate, the stud, and the mold. Furthermore, no heat trans- 
fer is assumed to occur across the boundary between the plate 
surface and the mold. Figure 8 illustrates this model 
configuration . 

The computerized model was developed to simulate 
temperature histories of points in the bar, plate, and weld 
metal as selected by a user. Since such information would 
be of long term use if underwater thermit welding were deve- 
loped into an operational process, the computer program was 
written to be as general as possible. In this respect, a 
user is able to specify different materials (but limited to 
those where exothermic reactions are possible) and different 
geometries, e.g., plate thickness and stud radius, in addi- 
tion to points for which temperature histories are desired. 

As stated in a previous section of this chapter, 
the model is based on finite difference techniques. The 
model starts simulation at the time the crucible is tapped, 
permitting molten metal to enter' the mold. The initial 
temperature of the molten metal is calculated based on the 
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INSULATION RADIUS 




Figure 8: MODEL CONFIGURATION AND FINITE DIFFERENCE GRID 
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theoretical reaction temperature of the particular thermit 
mixture and the device efficiency. Flow through the mold 
is provided for. 

Using these initial conditions and the environ- 
ment temperature specified by the user, the model calculates 
temperatures at all nodal points one time increment later. 
Then the model uses these latest temperatures to calculate 
the next future temperature, and so forth until the desired 
simulation time (specified by the user) is completed. The 
variation of material properties (conductivity and diffu- 
sivit.y) with temperature and heat release during solidifi- 
cation are accounted for in the computations. Temperature 
histories for the ten selected points are printed showing 
temperatures for each second of simulation time. 

Thus, the model steps its way through two indices: 
a space index and a time index. A simplified flow-chart is 
shown in Figure 9- The detailed-flow chart and the computer 

program, with sample input and output are included in 

* 

Appendix D. Appendix C presents instructions for program 



use . 
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Figure 9: SIMPLIFIED FLOW-CHART OF HEAT FLOW SIMULATION 
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III. EXPERIMENTAL PROCEDURE 

A. OVERVIEW: 

As discussed in Chapter I, the experimental 
approach is related to three objectives of this thesis: 
first, determination of the moisture level which can be 
tolerated between the surfaces to be joined; second, the 
degree of surface contamination which can be effectively 
removed by flow of the molten metal itself between the sur- 
faces to be joined; and third, verification of the model 
developed to predict temperature histories. In order to 
accomplish these objectives, a series of experiments was 
conducted utilizing a commercially available bar to plate 
thermit welding device. 

The material and equipment used, and the procedures 
followed are discussed in this chapter. 
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B. EQUIPMENT AND MATERIAL: 

The equipment and materials used are shown schem- 
atically in Figure 10 and are as follows: 

1. Thermit welding device 

2. Water tank 

3. Temperature recorders 
. Thermocouples 

5. Steel plate and bar 

6. Support and hold down apparatus 

A half-section of the thermit welding device used 
was shown in Figure 3 (Chapter I). The kit is designed for 
welding reinforcing bars to horizontal plates in air and 
is marketed by Thermex Metallurgical, Inc. of Lakehurst, 

New Jersey. This type of mold is fabricated using sand with 
a phenol-formaldehyde base resin as the bonding agent, and 
comes in two sections. The process is known as the "shell 
molding process." This type of mold is not reusable. 

Other materials included in the kit include mild 

* 

steel thermit welding compound (to match A.I.S.I. 1020 com- 
position), tapping disks, starting thermit, sealing compound, 
molding sand, and bolts (for joining the two half-molds). 

The thermit is reacted in a large chamber (on the 
left in Figure 3), and is tapped into the mold and between 
the surfaces to be joined. Flow of the molten metal con- 
tinues up a riser and into an overflow chamber (on the right 
in Figure 3) until the equilibuium level is reached or 
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EVENT 

MARKER 



Figure 10: EQUIPMENT SCHEMATIC 



solidification occurs. 



A steel water tank measuring four feet in length, 
two feet in width, and three feet deep was used for the 
experiments. The tank was- filled with ordinary tap water 
(Boston Ma . ) and drained by gravity discharge. 

Clevite Brush, Model 22, recording voltmeters were 
used to measure the output signal from the thermocouples. 
Three instruments were used for most experiments. Each 
instrument contained two channels so that a total of six 
thermocouples could be monitored simultaneously. Event 
markers on the instruments were connected electrically so 
that depressing a key on any one instrument would cause a 
simultaneous voltage signal to be recorded on all charts. 
Calibration was checked by inserting a thermocouple into 
known temperature baths (700°F and 1200°P). The instruments 
were adjusted so that chart velocity (horizontal scale) was 
1 mm/sec, and one division on the chart paper's vertical 
scale equaled one milivolt (or, 50 mv for total scale). 

Chromel-Alumel Type K, 0.020 inch diameter ther- 
mocouple wire was selected to measure temperature changes 
during welding. The size was considered more than adequate 
to withstand expected temperatures, and was governed primar- 
ily by the need for rapid response. In addition though, a 
small wire was considered desirable because thermocouples 
were to be installed by inserting and welding them into a 
hole drilled in the plate. Small holes were necessary to 



limit inaccuracies introduced into the problem. Installation 



is covered in the next section of this chapter. 

The plate and bar joined were mild steel (A.I.S.I. 
1020 composition) . The theoretical simulation model assumes 
a large plate. Therefore, the required length and width of 
the plate selected was based on the need to have little or 
no temperature change at the extremedies of the plate. A 
one foot square plate was expected to meet this requirement 
and was therefore selected. A thickness of one inch ( 4 0 
plate) was selected and considered satisfactory for verifi- 
cation of the model. 

The plate was supported by a platform suspended 
from the sides of the tank. The mold was secured to the 
plate (held down) by a horizontal bar resting on top of the 
mold and clamped to the sides of the tank. 

Ordinary mercury in glass thermometers were used 
to measure tank water temperature and the temperature of 
the thermocouple reference junction. 
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C. PROCEDURE: 

1 . Thermocouple Fabrication and Installation : 

In order to avoid temperature gradients, the 
thermocouple's hot .junction should be as small as possible. 
Ideally, it should be a point since we are, in effect, 
measuring a three dimensional distribution. Preparation 
of thermocouple hot junctions began with stripping two inches 
of insulation off the ends, and inserting them in 0.125 inch 
diameter, two hole ceramic insulators. The protruding ends 
were then twisted and fused together in an oxidizing oxya- 
cetylene flame until only a small spherical shaped (0.030 to 
0.040 inches in diameter) junction remained. Reference 
junctions were prepared in the same manner, only without 
the ceramic insulators. 

Since it was necessary to obtain temperature his- 
tories of interior points in the plate, 0.128 inch diameter 
holes were drilled from the underside of the plate at de- 
sired radii to the desired depths. The temperature measure- 

* 

ment locations are shown in Figure 11. 

Intimate contact between the metal plate and the 
thermocouple was necessary in order to avoid unnecessary 
measurement inaccuracies. This was accomplished by welding 
the thermocouple junction to the base of the hole in the 
plate. The welding method employed utilized a capacitor 
discharge apparatus. Leads were connected to the thermo- 
couple wires and the plate; the capacitor was charged; then 
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Figure 11: THERMOCOUPLE LOCATIONS 



the thermocouple was inserted into the plate to the base of 
the hole. The thermocouple was then withdrawn approximately 
1/16 inch, and the capicitor dischared. An arc was struck 
between the plate and the thermocouple junction (detected 
audibly), and the thermocouple was again inserted all the 
way to the bottom of the hole. In each case, a weld of 
sufficient quality u as achieved to obtain intimate contact 
between the thermocouple junction and the desired location 
in the plate. The joint was strengthened further by applying 
a ceramic cement (Sauereisen Sealing Cement) to the plate 
and ceramic insulator. 

Since the entire plate would be underwater, the 
cement was waterproofed with a five percent acetic acid 
solution, and a four foot length of shrink tubing was shrunk 
onto the thermocouple wire and ceramic insulator. Figure 
12 shows a section view of the completed thermocouple 
installation . 

% 

2 . The Welding Operation : 

A total of six experimental welds were performed. 

Of these, the first two were performed above water in a dry 
environment to ensure the equipment was functioning properly. 
The last four tests were performed underwater. It is the 
last four welds to which this section refers. 

Prior to each actual welding operation, several 
preparatory steps were necessary. First, the mold was 
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Figure 12: 



THERMOCOUPLE INSTALLATION 



assembled, and the tapping disk inserted. Since one objec- 
tive associated with the experiments was to determine allow- 
able wetness level between the bar and the plate, each mold 
was sealed against leaks so that water level in the mold 
could be controlled, and therefore, known. 

Second, the plate, with thermocouples installed, 
was placed on the supporting apparatus in the water tank. 

The thermocouple leads were then connected to the recording 
voltmeters (except in the last test). 

Third, the welding device was positioned on the 
plate and vertical pressure applied via the bar and clamp 
arrangement. The reaction chamber of the mold was then 
filled with four pounds of thermit mixture and a small amount 
of starting thermit. 

Fourth, the voltmeters were turned on and adjusted 
so that zero voltage was indicated while the plate was at 
ambient temperature (except in the last test). 

Fifth, the tank was filled with water until the 
plate was four inches beneath the water surface. 

Finally, after allowing sufficient time for the 
plate to attain thermal equilibrium with the water, the 
thermit was ignited with a oxyacetylene torch, and the 
welding operation begun. 

Temperature histories were obtained in all under- 
water tests except the last. In each of the first three 
tests, a different wetness level between the bar and the 



4 7 

plate was created. The first underwater test was conducted 
with no water in the mold; the second was conducted with the 
mold 25$ full of water; and the third was conducted with the 
mold 50$ full of water. In each of the first three under- 
water' tests, the surfaces to be joined were clean and free 
of oxidation. 

The fourth and last underwater test, for reasons 
which will become apparent in the following chapter, was 
somewhat unique. The experiment was conducted with the 
mold 50$ full of water, and on an oxidized plate surface. 
However, the major difference was the use of preheat. During 
device preparation, cavities were channeled in the device 
wall surrounding the mold. These cavities were filled with 
thermit which was ignited just before the thermit in the 
reaction chamber. The hypothesis was that by preheating the 
area to a high temperature, water within the mold would be 
vaporized and the actual welding would be performed in a 
nearly dry environment. No temperature measurements were 
conducted in this test. 

In order to ascertain overall quality of welds, 
the last three specimens were subjected to tensile tests. 



IV. RESULTS 



A. HEAT FLOW: 

In order to verify the computer simulation model, 
the points in the plate at which temperature histories were 
gathered matched those selected for simulation. Figures 
13a through 13f are plots of experimental data vs. theoreti- 
cal prediction curves. It is apparent from the figures that 
theoretical and experimental curves are quite close, thus 
validating the model as a prediction device. 

The greatest variance exists between theory and 
experiment for points "C" and "E." Note from Figure 11 that 
these points lie directly beneath the weld/insulation boun- 
dary. The simulation models this boundary as a perfect 
insulator. However, the experimental device did not offer 
perfect insulation, and heat losses through it existed. For 
this reason, theoretical simulation predicted higher temper- 
atures for points in close proximity to that boundary. 

The model could have been changed to account for 
heat losses through the insulation boundary, but it is 
hypothesized that an actual device would offer better insu- 
lation than the one used for test purposes. Therefore, 
simulated temperatures would be in closer agreement with 
those achieved with an operational device, and the model 
was left unchanged. 

Temperature history data tapes are included in 



Appendix E. 
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B. WETNESS LEVEL: 

Results with respect to an acceptable wetness level 
were somewhat disappointing. It was originally hypothesized 
that flow of molten metal through the mold would serve to 
evacuate any water present in the mold. This, however, did 
not occur. 

In review, the first underwater experiment was 
conducted with no moisture or water in the mold; the second 
was conducted with the mold 25% full of water; and the third 
was conducted with the mold 50% full of water. In each case, 
a weld, externally sound in appearance, was achieved. How- 
ever, when subjected to tensile tests, the second weld frac- 
tured at 25,000 pounds, and the third weld fractured at 
12,000 pounds. Interiors of each of these latter two showed 
extreme porosity. The first underwater weld (dry mold) 
showed no porosity. Figure 1*1 illustrates predicted strength 
vs. wetness level for a one inch stud to plate weld (mild 
steel) based on this data, and in the absence of any preheat. 

At 'this point I concluded that a sound underwater 
thermit bar to plate weld would not be possible unless a 
better method of evacuating any water in the mold were 
devised. The aspect of the experiment is discussed in the 
next section of this chapter. 
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FORCE (pounds) 




Figure 14; PREDICTED STRENGTH VS. WETNESS LEVEL 
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C. PREHEAT AND CLEANING ACTION: 

The final experiment was conducted using the same 
procedures used in other underwater welds except for three 
differences: temperature measurements were not taken; the 

area of the plate to be welded was oxidized; and preheat was 
used. As explained in Chapter III, cavities were channeled 
in the device wall surrounding the mold, one cavity in each 
half of the device. Each cavity contained sufficient space 
for 0.25 cubic inch of the thermit mixture. Holes were 
drilled from above the intended water line to the cavities 
and filled with starting thermit. 

After all other preparations had been made, in- 
cluding filling the mold half full of water, the thermit 
used for preheat was ignited immediately prior to igniting 
the thermit in the reaction chamber. Therefore, while the 
thermit in the reaction chamber was reacting, the water in 
the mold was being vaporized by the heat from the preheating 
thermit. Following the experiment, the weld was subjected 
to a tension test. At 44,400 pounds tension, the threaded 
portion of the bar fractured. This force was equivalent to 
subjecting the weld material to approximately 56,600 p.s.i. 

Rather than subject the weld to further destructive 
mechanical tests, the bar and plate were split, and a spe- 
cimen was polished and etched for visual examination. 

Figure 15 shows the entire weld area and portions of the 
plate and bar. The main features to note from the figure 
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are the lack of porosity in the weld and the excellent bonds 
achieved, even on the oxidized plate surface. 

Figures l6a, b, c, and d provide closer looks at 
the specimen at low magnification (50X). Figure 16a shows 
the weld material itself; Figure 16b shows the fusion line; 
Figure l6c is a photograph of the heat affected zone; and 
Figure l6d shows the unaffected base metal in the plate. 

Actually, there is evidence of a small porous area 
in the weld material. It is toward the upper right in 
Figure 15 and is shown under magnification (75X) in Figure 
17. The small amount of porosity in this weld, however, is 
certainly encouraging at this stage of research. 
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Figure 15: CROSS SECTION OF THERMIT WELD MADE UNDERWATER 

USING PREHEAT (PLATE AT BOTTOM, HX, 3 % NITAL ETCH) 
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Figure 16b: MICROPHOTOGRAPH OF FUSION ZONE 

(50X, 1 % NITAL ETCH) 
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Figure 16c : MICROFHOTOGRAPH OF HEAT AFFECTED ZONE 

(50X, 1 % NITAL ETCH) 




Figure l6d: MICROPHOTOGRAPH OF UNAFFECTED BASE METAL 

IN PLATE (50X, 1 % NITAL ETCH) 
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Figure 17: MICROPHOTOGRAPH OF POROSITY IN WELD METAL 

(75X, 3 % NITAL ETCH) 
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V. CONCEPTUAL DESIGN 
A. GENERAL DISCUSSION: 

As stated in Chapter I, the overall objective of 
this thesis is the conceptual design of a device which could 
be used to exothermicly weld a bar, stud, or padeye to 
underwater structures. The results presented in Chapter IV 
allow us to form important conclusions directly affecting 
device design. They are as follows: 

1. Water must be' removed from the weld area prior to 
tapping any molten metal into the area if a sound weld is 
to be achieved. It appears possible to accomplish this 
requirement through some means of preheating the area. 

2. Flow of molten metal through the mold can effec- 
tively remove minor oxidation and dirt from the surfaces to 
be joined. Since no data was obtained on the removal of 
paint or more serious contamination, it is assumed for design 

purposes that flow of molten metal through the mold will not 

* 

effectively remove anything more than minor oxidation and 
dirt . 



3. Finally, and possibly most important, the success 
of the experiments demonstrates the feasibility of the pro- 
cess, thus implying that time spent on conceptual design and 
further development is not fruitless. 



Device design is also guided by some basic 
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assumptions. The Navy's interest lies in salvaging or 
raising an object from the ocean floor. The device, there- 
fore, should be operable at different (and deep) depths. 
Since the depths may exceed those reachable by a diver, the 
device should be remotely emplacable and operable through 
the use of a submersible. Furthermore, it is unrealistic 
to assume that the object has a horizontal surface on which 
to weld. For this reason, any device should be capable of 
operating at angles from zero to ninety degrees. 
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B. PRELIMINARY CALCULATIONS: ' 

1 . Size and Strength : 

The original specifications of the Navy were that 

the padeye have a load carrying capability (L ) of 50 tons 

Y 

or 112,000 pounds. Using a thermit mixture of mild steel 

composition assuming a yield strength (a ) of 60,000 p.s.i., 

«y 

the following calculations for the cross-sectional area (A) 
of the weld can be made. 

L y ■ A ’°y 

A = L /o = 112,000/60,000 
A = 1.868 in. 

If the area is circular, this equates to a radius of 0.772 
inches. Based on these figures, the load carrying portion 
of the total configuration was designed and is illustrated 
in Figure 18. The actual diameter of the area to be welded 
is 1.75 inches. Based on weld material strength of 60,000 

p.s.i., this configuration would have a load carrying 

* 

capability of 131,000 pounds or 58.6 tons. 

2 . Preheat : 

A rather "make-shift" means of providing preheat 
was employed in the experimentation. Although special 
thermit mixtures exist specifically for providing preheat, 
previous research has indicated that another method may 
serve to not only accomplish this purpose but also remove 
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Figure 18 : LOAD BEARING PORTION OF THERMIT WELDER 
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contaminants such as paint and minor marine growth from the 
plate surface. 

This other method utilizes rocket technology.^ 

The basic concept is that by exhausting high temperature 
gas into the weld area, both water removal and cleaning can 
be accomplished. An additional benefit is that some pre- 
heating of the parts to be joined is also accomplished. 

The design is based on existing propellant compo- 
sition, manufacturing methods, and assumptions as follows: 

a. Propellant: Arcite 4l3A 

b . Wired end burner 

O 

c. Mass flow rate (M) = 1.7 lb/sec 

d. Density ( p) = 0.061 lb/in 3 

e. Chamber pressure (P ) = 2500 psi 

f. Burn rate (wired) (r^) =1.8 in/sec 

g. Gas temperature (Tg as ) = 5000°F 

h. Ambient temperature ^ ) = 40°F 

i. , Desired burn time (t^) = 3.0 sec 



Given these specifications, grain diameter (D ) and grain 

S 

length (L ) can be calculated as follows: 

o 



Q 

See Masubuchi, Koichi, Final Report on Underwater 
Thermit Welding of 50-Ton Padeye , Massachusetts Institute of 
Technology, October 29, 1971. 

•^ Ibid . , Appendix C, p. 4. 
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D 



g TT /4) *P* r b 




M 



O 




IT • 0 . 061*1. B 



4 • 1 . 7 



D 



= 4.42 in 
g 




L 



g 



5.4 in 



A grain of these dimensions weighs 5-05 pounds. According 
to Atlantic Research Corp., such a grain is easily manufac- 



Present thermit welding devices used on the surface 
rely on gravity induced flow of the molten metal from the 
reaction chamber into the mold. Since in a device designed 
for underwater use the thermit reaction takes place in a 
chamber isolated from the environment, the molten metal must 
be forced into the mold area. None of the theoretical pro- 
ducts of the thermit reaction are gaseous, and therefore, 
the pressure developed in the chamber upon reaction depends 
upon the temperature, pressure and volume of gas which was 
in the chamber prior to initiating the reaction. Although 
no specific calculations are presented here, the reaction 
chamber pressure requirements can be fulfilled prior to de- 
vice emplacement based on temperature, volume and the speci- 
fic- gas used . 



tured within the present state-of-the-art. 
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3 . Pressure : 



11 Ib_id L . , Appendix C, p. 4. 
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C. DESIGN: 

A conceptual design was prepared based on results 
presented thus far in the thesis and is illustrated in 
Figure 19. The upper portion of the figure illustrates the 
device as it would be viewed from above while resting on a 
horizontal plate. Only the exterior outlines of the rocket 
propellant and thermit reaction chambers (right half of the 
view) are shown in this view. The lower portion of the 
figure shows a side section view. 

The major portion of the device could be con- 
structed from mild steel. The reaction chambers are of 
equal size and inclined at an angle of 45°. Each chamber is 
lined with a refractory such as magnesia. Ignition timing 
control and thermit chamber pressurizing apparatuses are 
located between the two chambers. Void volumes of the de- 
vice are filled with insulating material and ballast to 
achieve the desired weight in water. 

Perhaps the best way of explaining the device 
would be a discussion of its sequential operations. Before 
emplacing the device, and while still on the surface, the 
rocket propellant and thermit (6 pounds) are placed in the 
proper chambers, and the chambers sealed. The thermit 
chamber is then pressurized to a level dependent upon the 
depth of the intended weld, so that when the thermit reaction 
is complete, pressure equilibrium between the chamber and 



the environment will exist. 
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Figure 19: UNDERWATER THERMIT WELDER 
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The thermit reaction must be initiated first. 

Just prior to completion of the thermit reaction, the rocket 
propellant should be ignited. This timing and overall delay 
are controlled by the ignition timing control apparatus. 
However, at this writing, times reported in the literature 
to complete the thermit reaction are not accurate enough to 
state specific sequencing times. See Chapter VI for research 
recommendations. Suffice it to say at this point that the 
ignition timing control would be set and sealed. 

The device is then emplaced by a diver or submer- 
sible with the reaction chambers on the "uphill" side if 
the object's surface is at an angle other than horizontal. 
Special appendages for handling purposes are not shown in 
Figure 19 , but if attached, would not interfere with device 
operation. Temporary means of attaching an operational 
device until the weld is completed were not investigated and 
should be considered in any further study. The strength of 
any temporary bond need not be great, but only sufficient to 
hold the device in place until the reactions and the weld 
are completed. Adhesives or small thermit tack welds may 
suffice i 

The timing mechanism first ignites the thermit. 
Fifteen to twenty seconds later, the rocket propellant is 
ignited, evacuating the water in the mold and vent riser, 
and preheating the area. When the propellant chamber pres- 
sure is in equilibrium with the environment pressure, the 
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molten metal Is tapped Into the mold. A portion of the 
metal continues up the vent riser and solidifies as it con- 
tacts the water, thus isolating the mold area from the water. 
The remaining metal seeks equilibuium between the propellant 
chamber and the thermit reaction chamber. 

The gating is sized so that molten metal in the 
mold is always flowing under positive ferrostatic pressure. 
Since the chamber axes are at 45° to the object surface, 
this pressure will exist when welding at any angle from 
zero to ninety degrees. 

Less than five minutes after the thermit reaction 
was initiated, the device should be capable of bearing the 
prescribed load. Due to time constraints, a test device 
could not be constructed, developed and tested as part of 
the thesis. Recommendations are contained in the following 



chapter . 
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VI. CONCLUSIONS AND RECOMMENDATIONS 
A. CONCLUSIONS: 

Perhaps the most signicicant conclusion arising 
out of this theses is that sound welds can be produced 
underwater by utilizing the thermit process. This was dem- 
onstrated by successfully welding a bar to a plate under- 
water. - ■ ‘ ' 

However, several otner important conclusions can 
be made. First, the accuracy of the temperature simulation 
model would seem to validate the hypothesis about heat losses 
underwater. More specifically, the boiling mechanisms on 
the bottom of a horizontal plate do not provide nearly as 
much heat transfer as boiling on top of the plate. Radia- 
tion through the vapor film would be the primary mechanism. 

Second, a means of removing the water from between 
the surfaces to be joined must be incorporated in any thermit 
welding device destined for underwater use. Unless this is 
accomplished, porous welds can be expected. 

Third, flow of molten metal through the mold can 
effectively remove minor oxidation and dirt from the surfaces 
to be joined. 
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B. RECOMMENDATIONS: 

The first recommendation is related to underwater 
welding in general. If temperature simulation models are to 
become useful devices, more accurate definition of heat loss 
terms is essential. It is recommended that basic research 
be initiated on convection and pool boiling generated by a 
localized heat source on flat plates at all orientations. 

The other recommendations contained herein are 
related specifically to underwater thermit welding. A 
requirement is envisioned for both academic and develop- 
mental work. Recommendations are as follows: 

1. Construct and develop a device similar to that 
discussed in Chapter V. Several associated problem areas 
can be identified. 

a. Accurately define speeds of thermit reactions 

b. Develop timing and ignition devices accor- 
dingly . 

c . " Confirm sizing and arrangement of gating. 

d. Investigate methods of temporarily securing 
the device to a surface. 

e. Investigate dynamics of molten metal flow 
and pressure balances. 

2. I believe the thermit process to be feasible for 
other underwater applications such as joining structures, 



72 

cables and pipe lines. Given this assumption, long term 
effort is implied. Development of other devices is essen- 
tial. However, a most important research area lies in 
improving the properties of underwater thermit welds so 
that specifications of appropriate organizations and 
governing bodies can be met. 
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APPENDIX A 



PARTICULAR SOLUTION TO HEAT FLOW EQUATION 

The general equation is: 

V 2 T = I.H 
a at 

Assuming temperature can be expressed as a product of time 
and space variables as follows: 

T = U ( r , s , 0 ) • 4> ( t ) 

Then : 



V 2 T = V 2 ( U * 4> ) = cj) • V Z U 

and 

9T = Ti.ii 
3t u 8t 



Therefore , 

• v 2 u = I.u.|f 

v a ot 

and 

III = _I_.M = _ Y 2 

U a • <j> 3t Y 
where y 2 is a constant. 
Since $ is a function of t alone, 

3d> _ d<b 
3t dt 

and 



44 + y 2 'a • <j)= 0 
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to which the solution is 



<t» = A e 



-Y at 



where A is a constant 



Also , 



V 2 U + y 2 U = 0 



Assuming symmetry with respect to 0 in cylindrical coordinates 



V 2 U + v 2 U = 



3 2 U + I.3U + 3 U 
3 r ; r * 3 r 3 z 2 



+ y -U = 0 



Again assuming a product solution, 



U = R(r ) • Z ( z ) 



Since R is a function of r alone, and Z a function of z 
alone, it follows that: 



3 2 R , Z . 3R 7 „ . 3 2 Z 9 .,2 



+ — • — — + R • 



3r 2 r 3 r 



3 z 2 



+ y *R*Z = 0 



and 



1 d 2 R , 1 dR ,1 d 2 Z 



R dr 2 R • r dr 

-r 2 : - - 



= _(±._ — _ + y 2 ) = 

' r ? rl r? 2 Y . / -> 



Z dz 2 
constant . 



Therefore : 



■ 2 — y + r + £ 2 *r 2 *R = 0 
dr 2 dr ^ 



to which the solution is; (for 5 / 0 ): 



R = ZoUr) = B • J 0 ( £r ) + C'YpUr) 
where B and C are constants. 
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and , 

0 + (Y 2 - C 2 )‘Z o o 

to which the solution is: 

Z = D‘e ~ iz>/ ^ 2 - E-e lz/lf2 _~ ^ 
where D and E are constants. 

Thus, the particular solution for temperature (T) is the 
product of the separated solutions: 

T = <fr(t ) -R(r) • Z ( z ) 

or, . 

T = [A • e~ Y at ] -[B-J q U r) + C«Y 0 Ur)] 
rn -iz/y 2 - C 2 , p iz/y 2 - c , 2 

where A, B, C, E, E, y and £ are con- 
stants which must be evaluated by 
applying the applicable boundary and 



initial conditions. 
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APPENDIX B 

DERIVATION OF HEAT LOSS TERMS 

1 . Top of Plate Convection: 

The following equation is recommended by McAdams 
for convection when the top surface of a horizontal plate 

• v, +- 12 

is hot: 



~t~ ' O.HWG^-PrpVa 

All properties for convection purposes were taken at 75°F 
and are as follows: 

k £ = 0.35 [BTU/hr/ f t/°P] 

Pr^ = 6.4 
Gr 2 

— T ■■■■■■ = g * s * - p . - = 4 . 4 • 10 7 [l/°P/ft 3 ] 

L -AT u 2 

where: g = acc. of gravity 

P = density 

, 6 = coef. of expansion 

y = viscosity 

L = length 

AT = T - T 

plate water 



"^See Kreith, ojd. cit . , p 340. 
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Substitution of these values into the convection equation 
yields : 



or 



h 



c 



32.2* (T 



plate 



T i )V 3 
water 



9 C /A - 32.2-(T Dlate - T Mater )V 3 [BTU/hr/ft 2 ] 



2 . Bottom of Plate Convection: 

The recommended equation in this case is:^ 



h • L 
c 

k„ 



All properties 
result : 



= 0 . 27 * (Gr^-Pr^) V 4 

are again taken at 75°F yielding the following 



h = 12 2 • (T -T ) 1 / h 

c v plate water' 



or 



q c /fl = 12.2-(T nlate - T water )V» [BTU/hr/ft 2 ] 



3 . Top of Plate Nucleate Boiling: 

The correlation presented by Rohsenow and Choi 
is as follows: 

V &T . , V A * / g °' q 1 

h f g - pr f 7 " sf *v h fg 

Pq See Kreith, 0£. cit . , p. 390. 
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Solving for q^/A: 



V A = 



14 



V AT 



C •T) • Pr 1 • 7 
v °sf n fg ir £ 



u r h f K 



go -a 



g-(Pj,-P v ) 



Properties for nucleate boiling are taken at the saturation 
point (14.7 psia, 212°P) and are as follows: 

C £ = 1.021 BTU/lbm-°F 
c sf = 0.013 
hfg = 970.3 BTU/lbm 
Pr £ = 1.9 

g 0 = 4 . 17 • 10 8 lbm-ft/lbf-hr 2 
g = 4 . 17 • 10 8 ft/hr 2 
p £ = 59-8 lb/ft 3 
o v = 0.0372 lbm/ft 3 
a = 0.00404 lbf/ft 
y £ = 0.74 lbm/hr-f t 



Substitution, of these values yields the following result: 
%/A * 1 -7‘t9MT Dlate - T sat U [BTU/hr/ft 2 ] 



11 ) 



See Rohsenow and Choi, op. cit . , p. 225. 
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APPENDIX C 

INSTRUCTIONS FOR PROGRAM USE 

The temperature simulation program has been written 
in FORTRAN language and can be executed without difficulty 
on any computer having a FORTRAN Compiler and 153>600 bytes 
of primary storage for use by the program. The user is 
required to furnish the following information describing the 
material, the environment, the device configuration, and the 
points for which temperature histories are desired: 

1. First data card ( FORMAT ( 4F10 . 4 )) : 

TENV : temperature of the water where the weld is 

to be made. [°F] 

TSAT: water saturation temperature at the depth 

the weld is to be made. [°F] 

TSOL: solidification temperature of the metal. 

[°F] 

THERM: theoretical reaction temperature. [°F] 

* 

2. Second data card (FORMAT (5F10.4)) 

GAP: distance between the bar and the plate. [in] 

REAR: radius of the bar. [in] 

RINS: radius of the insulation. [in] 

THICK: plate thickness. [in] 

EFF: theoretical or experimental device efficiency. 
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3. Third data card (FORMAT (2F10.il) ) 

RHO : density of metal at fusion temperature. 

[lbm/in 3 ] 

QFUS: latent heat of fusion of metal. [Btu/lbm] 

4. Fourth data card (FORMAT (13)) 

IVAL: the number of data cards Immediately 

following which describe thermal properties of 
the plate and bar. (maximum of 100) 

5. Fifth group of data cards (FORMAT (3F10.5)) 

TPROP: temperature at which particular property 

holds. [°R] 

CON: thermal concuctivity at temperature. 

[Btu/ft-hr-°R] 

DIF: thermal diffusivity at temperature. [ft/hr] 

6. Next data card (actual number depends upon the 
number of cards describing properties) (FORMAT 
(2F10.4)) 

TSIM: desired simulation time. [sec] 

FLOW: time of flow through the mold. [sec] 

7. Next to last data card (FORMAT (1215)) 

Cylindrical coordinates of six points in the plate 
for which temperature histories are desired in the 
number of eighths of an inch a point is away from 
the center of the weld and plate top. 
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8. Last data card (FORMAT (415)) 

Coordinates of four points in the bar or weld 
material for which temperature histories are 
desired in the number of eighths of an inch a 
point is away fron the top of the plate. 



APPENDIX D 



PROGRAM FLOWCHART 
PROGRAM LISTING 



SAMPLE INPUT AND OUTPUT 
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0 * /? 7 / 7 2 
CHART TlTlf 



AUT jFI OW CHART SC T - T» “ ?r •» ATH» «• . Jm j L *f t«lN )GP PA*, 01 



intpcouctjry cements 



PROGRAM PCR PRC01CTJNG TE“Pr° ATOP E OT T I s' » 4A u T 
PI AT ^ THERMIT WFIO UMn"PWA T rB 



(HAH Sf T - 



• MOL AT | IN Pf ,P 



i n 



/ VI A F t- MAT 
f I *m*% / 

/IN* ’M t I St / 



/ v: 
/ :\t 



I os 

/ P r Af fH M ncv / 
/ 5 / 

/ Vlt F'iPMAT / 

/ IOOP / 

/ INTC THE l 1ST / 

I 

I 

I NOT-* Of 
*■ LT ST = FHT, qFijs . 



/ HCAf ff~m ntv / 
/ e / 

/ yli fifiMiT / 

/ IOOh / 

/ INTC Thf I I<T / 

I 

I 

I not; oh 
* list * ival • 



/ VIA r B*AT / 
/ IOO* / 

/ t n r " r h‘ i r/ / 

l 

I 

I N. T 

* LIST a 

« I fPi '0(|Ptf“'P|, 

* r m ip-' p>. 

* 0IM1 w < io| « 

* t pa p =• l « IVAL I 



/ IM I F P f'M OFV / 

/ c / 

/ vl ' r umat / 

/ 100" / 

/ l'.T • L J c T / 



I NCT 12 



21 

l w* 0. r I 



IF V 
IHXI C V 



IFtxML 

THi-PM = 



f /Oe L r I 

M - 

w/rULT I 
Ttl'SM - 



>• 

I 





/ Ft 1 A(J Ft- M 


O' V / 


l 

1 




1 OFtJS » - CPUS 


1 




/ S 


/ 


1 




• 






/ VI A F JMM4T 


/ 


1 




i 




/ 


ton 


/ 


1 




i 




/ 


1 NT THt l | ,r 


/ 


1 




1 NOTE 


z* 








1 




• •••••••• 


• * 




I 




1 




• BEGIN 00 LO »p 


• 




1 




1 




• 25 1/ = I, 10 


• 




1 N' 7 >- 1 A 


1 








- 


t ! r t r 1(4, 


JA, • 


1 




1 






! t j-. r r . 


r * * 


1 


I 


1 


25 


• 


! ?' . J . I I , 


j . * 


1 


1 






* 


1 I F Jf 




1 


1 


1 BET *(121 - 999.9 


J 




J 




1 

1 


1 

1 


1 LTJlin i 0.0 


1 




/ t Tv 


vv / 


1 

1 


i 


I 

1 






/ 


/ 


1 


i 


25 * 26 






/ V . 4 - « •• T 


/ 


1 


i 


• • 




/ 


10 l l 


/ 


1 


i 


• • 




/ 


• ** •«. | 


/ 


1 


i 


NO • • 










1 




• ENU JF O'* 






1 




1 




• L JOP? * 






1 




1 




* • 






1 •)< 


T* If. 


1 




• • 










1 








• 


LISI * 1 r, . I h , • 


1 




IVES 






U. M 




1 

1 




1 

1 






1 

1 


17 


1 

I 




1 

1 










1 




1 


27 




/ 


/ 


1 










/ WRITE TO OEV / 


1 




1 TENO TSIM/OELT 


1 




/ * 


/ 


1 




1 


1 


/ 


via format 


/ 


1 




1 LTEXiO * 


1 


/ 


1000 


/ 


1 




I IF IX 1 TENO ) 


1 




l 




1 




1 ITDP =• 100 


1 




1 


19 


1 




1 ITOPB * 99 


1 




/ 


/ 


1 




1 






/ WRITE TO OEV / 


1 




1 






/ A 


/ 


1 




1 


2? 


/ 


VIA FORMAT 


/ 


1 








/ 


tool 


/ 


1 




1 1 1 ROT = 


1 









1 




1 IF I XITHICA/SJ ♦ 1 


L 1 




1 




1 




1 JMAX • A 9 


1 




1 


i<* 


I 




1 


1 







• 


1 




| JMAXP r SO 


1 


1 


PI = 3. 1 A 1 5 9 I 


1 




1 


1 


1 




1 


1 




1 . J8AR * 


1 


1 


S = 0.125 


1 


1 

I 




1 IFIX(PBAR/SI ♦ l 


1 



1 

I 

I 

JBARM * JBAP - l 
JINS = 

IF I X I P I NS/S I ♦ l 

l PP * 

IF MCI .0/0 LTl 



* R C GI M O'J LOOP 

* 10 J Hi JHAXP 

I 

| < 

I 

I R(J J = Rl J - 1) ♦ 

I S 



END OF 00 
LOOP? 



ITE 5 

I 



BFGIN 00 L OOP 
20 I ■= l, IT C® 



| K 3 a 



UUI 
UP I I I 



TENV 
C eNV 



E NO OF 0? 
LOOP? 



nrcs 

I 

l 

I 

I 

I 

I NOTE 

• BEGIN 00 LOCP 

* 30 II = l, IIBoT 



02. AO >1 

I N ?T E 

• BEGIN O'J LOOP 

• 30 i - l. J“AXP 

I 

l< 

I 

I T1II.JI ^ MSV 
I 

I TP( I I , J) = TEKV 



I YE S 

I 
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CM! / ? 



*nr fi’w t hast sr t 






PAG. 0» 



cufioi i|T| _ <>p - r tout f $ 



02.40 — >• 

I 31 



| [UP j |T1P - | 

I I F f M f T.A P/ S I I 



I 

I 

I N'.TE 02 



• Pf0 I N 03 LC IP • 

• <-0 I • IGAP, 

• hop* • 



I 

>1 

I 03 



e NO ’F 01 
l )P ? 



1 1 A * ! I A • 
JA • JA * 
IIP * HR ♦ 






IIC 



JO 



IIC ♦ 

1 

I 

I 

JC = J" ♦ I 

no = no ♦ 



JO 
I IE 
JF 



jn 



i 

l 

I o? 

I I I F = 1 1 F * I I 

I I 

I JF * JF * 1 | 

I I 

I IG = 100 - IG 1 

I I 



I IH * 100 - Ih 




I 

I 01 



I I J = 100 - I J I 

I I 

I IK = 100 - !« I 



I 

I 

1 

ST AO T T I t*F L.'*lP 

I 

I 

I N ITf 10 



• PEG IN 00 irr.P * 

• TOO t * 1. I TEND • 



I 

OF. 21 > I 

I 

PAP T ( **P C At CULAT ION 

P'lL 

I 

I 

I 

I / 



/ 4.01 
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04/ 2/ 



AIITPFL w C HAP T SET 



rt**pf DA FU-> f- SI*lHAT|n«* PRjf,* 



P Af,F 04 



r HAM M.T l f - fUK tniJP F* 



C 3. IC >♦ 

I V ?«• 01 



* I* . I *1 r l • L I ' o 

♦ 100 I * *,0. I T 'pH • 



C c .Q-7-- 



I IY = I - JG1P 

I 

I 

-»0l I 



lire ip - CUII • I 

1 1 ♦ i H ♦ n ♦ i 

bllll/UO ♦ 4*0.0 I 



I PF P =. 

IF I * C IIP**' o/« 0 . 0 ) 



1 PPJP .GT. 
* ! VAL 



ITPUE 

! 

I 

I 

I 

• 

I 

I WOP * l VAL 



I 

->l 

I 



DIFUS * 

OIF( I PPCPI /2S.0 

CC?NDU = 

COM I POO*) 

/4 3200 • C 



S*»2/tOTFUS*OELTI 

I 

I ( 



I IJPII) - lUlf - 
I II * on ♦ 

I 111/* ♦ 11.0“ 

1 2.0/M|*Utn 



T .IT. ICAP 



FALSE . S 

. 0 






l . GT . I Fl W 



BEGIN 
99 lw 

in 



I LOOP 

IGAP, 



END OF PO 
LOOP? 



I YES 

I 



• L .GT. LC VA r 



STA3 T sn L 101 F ICATlfN 
R’TIJT INt 

> l 

997 • 17 



• U(I) .GE 
TSOL • AND 
•UPI1I ,L E 
* TSOL 



IFALSE 

I 

I 

I 

I 

I 

I 



>• 

I 



PETAIIYI . 
• <>99.9 



• I 

• • • • • 
IFALSE . S . 
I . 01 . 

I .... 

| 996 

I 
I 
I 

• 19 



BfcTAUYI • TRUE 
.LT. 11.0 - *-♦ 

ETAIIYII • 1 

• • I 

• • I 

IFALSE . S . 
I . OS . 

I .... 

I 994 

I 



I 



| Nl T = 

• eiGJN D} LOFIP 

• I* * I GAP, 

• irjp* 

I 

I < 

I 

I UPIIXI * Them 
I 
I 



FNO Jf 00 
LOOP7 



I YES 

I 

I 



too 
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(V* H ' / 7 2 



PUT 'H 1 w LMA4T 



T( HOI k A TUP I » (Mill AM'S P« ,r,k 



PA<> 



“A 7 Tt ft; - PU >1 t >u«f ^ 



/ / 






/ >54 / 




1 

On. IP >| 






1 

04 . p — >i 




1 


01 




i 


04 


1 i|p< it . T S' ‘L 


1 




1 ETAtlYI - 


1 


1 


1 




1 P T At | Y J » 


1 


1 £Y A < l r ( - (IJ( | . 


1 




1 RfTAl I Y 1 


I 


1 11 ♦ till ♦ 11 - 


1 




• . 




l ?.o* nit - 


1 




I 




I M * t T S >L - 










1 Ul I 111 


1 1 




I 


04 


1 Ang.-^ T / ( : h , 


t 1 








1 CFUS*<-**?I 


1 1 




1 8 M»(IY 1 * tun - 


■ 1 








l n * ut t ♦ i » - 


1 


1 


1 




i 2 . 0 -itm 


I 


I 


t 




1 • r*LT*' NOU/ 1 RMl* 


1 


1 


1 




1 JPljW*?) 


1 


• 0 ? 


1 








• • 


1 




I 




• • 


1 




-o-. ;••—>! 




* • 


Tktl | 


1 


t 




* f T A* lv » .L T . 


» * 


1 


lO • 07 




• 1.0 * 




I 


• • 




• * 




1 


• . 




• • 




1 


• « 


N r 


• 




1 


* no <r • 


_ » 


If ALS 
t 
| 




1 

1 


• L p? • 


1 

1 


1 

1 




1 

1 


1 Y 




1 

1 




1 

| 


1 


02 . 


• 03 




1 


r 




• • 




1 


i 




• • 




| 


i 




• • 


T RUE 


1 


i 




• FTAt I r t ,E 0 . • 





A 


t NO eAP TE NPE PATURE 




• 1 .0 * 






CALC UL ATI JV 










1 




• • 






1 




• 






1 NOTE 


08 


IFALSE 










1 






• BEGIN 00 LOOP 


• 


1 

1 






• 200 I t = 1 , l 133 T 


• 


1 

Oh. L*~-> | 






1 

08 .07 >| 




555 1 


04 




1 NOTE 


0 ° 


• 


- • 




• • ••««*••• 


• 


» uptn = luti - 


1 




• 8 e G IN 00 L 00 ° 


• 


» 11 • U I 1 * 11 - 


1 




• 2 CO J = 1 . J *AX 


• 


| I P Hi* 3 PuS *S* * 2 1 


1 








1 /C 3 N 0 UI / P ♦ 


1 




1 




1 < 1.0 - 


1 




08 . 06 — -> | 




1 i.O/HI*llIII 


1 




1 






1 




• 10 




1 ETA I I Y t * 0.0 


1 




• • 










• • 




1 






• 11 .' 0 . I • 


TSUc 


1 






• .ANO. J .LE. • 


— 



1 • J 8 AR • 



. 5.07. 

... 100 



I 

1 

I 

I 

I 

1 



US? 



A V" T£HP AND |0 ! N ▼ » 
AVc CALCULATION 
POLL 0 *S 



*>0 I II 



AREA * | 

P 1*1 S/ 2.0 \ 

I 

TAVE - I 

TI 2 ,n«ARPA | 



I 

I 

I NOTE 12 



• BEGIN 00 LOOP • 

* 300 K = 2 , J 0 ARM • 



I 

>| 

1 13 



A PE A ■ ARE* ♦ | 

2 . 0 *P ? *R ( A I *S I 

1 

T AV C = T A V - ♦ I 
T« 2 .k» I 

• 2 .o*Pi*ffim*s i 



1 

I 

I 

300 • 14 



NO • • 

ENO OF OC • 
• LOOP -7 • 



I YES 

t 

I 

I 

1 

I 

I 15 



AP C AR = 
PIMPIKt - 
S/ 4 . 0 )*S 

TAVc * TAVE ♦ 
T I 2 , J 8 AR ) • APR AR 

AREA = AREA ♦ 
AREAP 



I 

I 

I 16 

T AV E * TAVE/AREA | 

LlPPfjP = I 

IUI 1 T 0 PM) ♦ | 

TAVE • 1 

Tt II, m/ 3.0 ♦ I 

<.{>0.0 I 

TPROP = I 

lFlXtiJPP^P/ 50.01 I 



I 

1 

I 

• 17 



• • FAL ! 
• 1 PR 0 P ,GT. • 

• I VAL *1 

* * I 

. 1 TRUE . 5 

I . 1 « 

I 

t 

1 

I 

I 18 

I IPPJP • I VAL I 



05.17* — >• 

I 15 

- — — • 

OtFUS = I 

OIF I IPROP I /2 5.0 I 



S**2/tOIFUS*DELT» 

TP( I I ,JI = 
tUItToPNl ♦ 
TAVE1/1 ♦ tl.O - 
2.0/* t *T( 1 I , J) 



I 

I 

I ^0 

• 

I IJP t I T r P ) = I 

I TPH.ll I 



I 

I 
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8 . 06 , 



200 



04 tat }/ 



ajTrfj. 4 t m a p t sft 



Tfxfrt oafuo- s IHULAT |ik ph (,p 



t > b r .t oo 



I*H4 vl T» Tl 4 - OK I fniiK 



/ 9 I 



£'' M« )l iT| N 

c l l 

C*>. !<)--->♦ 

I 

«0 • 01 



I 

->l 

I 



04 



I or *. nv * - 

I .00306® l • 4BS I T 
I (II, tl - 

I Tt NV ( 4 * 1. 33033 

I OB'ML * 0.0 





I 

I 

I 



OB 



J .r,r. jtr 4 s 



TauE • • 

* TIII.JI .LE. • 

i • tsar 



|F»lSF 

I 



I 

t 0*> 

I OPftl = - | 

1 .O000014««eSIT | 

1 III,))- I 

I TS«T|*,1.0 I 

I 

>1 

PI I I 07 



I QTPTl - QBOIL ♦ I 
I •OCk, NV | 



I 

>1 

32 I 09 

uch:p * mu ♦ I 

I.J)*R(J) ♦ I 
TI1WJ - 1 1 • R I J - I 

n ♦ Tiii.j ♦ i 

I I • R ( J ♦ 1) ♦ I 
TUI,J|*R(JI> | 
/( 2 • 0 »R ( J I ♦ I 
R(J - II ♦ P l J ♦ I 
11)4 <>*0.0 I 



I 

I 

I OR 



| IPK-P * I 

I IF! X( UPRTP/50. 01 I 

I ♦ I I 



I 

I 

I 



10 



f 41 SE ■ 

• iprop .GT. * 

) • IVM. • 



I TRUE 



I 

I 




I { PR HP = TV4L I 



I 

>1 

I 12 



OI*=<JS * 

0 1 F I IPR3PI/2B.O 

C r NOU * 

C M I P p "'P I 
/ 4 3 20 0. 0 



S**2/(0I p US*n£L Tl 



I 




>• 

I I 3 

TP ( I I , J ) a I 

<2.o*n n ♦ i 
i ,ji ♦ tiii.j - i 
i * •« i.o - l 
s/i2.o*mj)ii ♦ I 

fiitj ♦ I 

Il'l 1.0 ♦ I 

B/I 2 . 0 -PC Jl II * I 
IQT_TL/ r ' l NDUI I 

•2 .0* S I ♦ I 
I 1.0 - I 

4.0/H|4 THI ,JI I 



I 



I 



89 



b.O*. 



0 ■ / l U * i 

r HAII T ) t (. F 



AIJT >n M f H A P T St T 



7 F^RMATUi- ► 5 I**UL A f I JN pp jGP PA', 



2u t UTl H|-p c LI 

Oj.n t — >• 

I 

•o * 01 



20 IN B J T TC 
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r i ."o. i I’n 



I 

I 

I 

I o 

U RR*1P • mi! - 
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